Steering Bose-Einstein condensates despite time symmetry 
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A Bose-Einstein condensate in an oscillating spatially asymmetric potential is shown to exhibit 
a directed current for unbiased initial conditions despite time symmetry. This phenomenon occurs 
only if the interaction between atoms, treated in mean-field approximation, exceeds a critical value. 
Our findings can be described with a three-mode model (TMM). These TMM results corroborate 
well with a many-body study over a time scale which increases with increasing atom number. The 
duration of this time scale probes the validity of the used mean-field approximation. 
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The realization of Bose-Einstein condensates 
(BECs) [l| of dilute atomic gases has opened new 
possibilities for the investigation of interesting physical 
phenomena induced by atom-atom interactions [3]. For 
example, it was found that for a BEG in a double- 
well potential the tunneling can be suppressed due to 
interaction- induced self-trapping 0, 0, Q, thus main- 
taining population inbalance between the two wells. 
A similar phenomenon has been predicted also for a 
periodically driven BEG in a double well @, 0, S, H, 
and shows the importance of interaction in the coherent 
control of quantum tunneling between wells. 

In recent years, ample interest has arisen in the field 
of directed transport in classical and quantum systems 
in absence of a dc-bias, see 11, Il31 and Refs. therein. 



In particular it has been shown that by breaking space- 
and time- inversion symmetries it is possible to achieve di- 
rected (ratchet) transport [l3, 14|. The study of the role 
of interactions on the quantum ratchet transport, how- 
ever, truly remains terra incognita, with only a very few 



preliminary studies available [H, 3, 12 1- For instance, 
in [isl it has been shown that atom-atom interactions in a 
BEG can cause a finite current while the non-interacting 
system would not exhibit directed transport; however, all 
these systems possess a Hamiltonian that is both space- 
and time- asymmetric. 

With this study, instead, we consider a non-dissipativc, 
i.e. Hamiltonian system in which the time-reversal sym- 
metry is not broken and the system is driven smoothly. 
We show that, starting out from initial conditions 
symmetric both in momentum and space, interactions, 
treated in the limits of a mean-field approximation, in- 
duce directed current in a BEG. In a full many-body anal- 
ysis of the system, symmetry considerations imply that 
this directed current is asymptotically decaying. Despite 
this fact we show that a finite current persists over a time 



scale which increases with increasing atom number in the 
condensate. 

To start out, we consider a BEG composed of N atoms 
confined in a toroidal trap of radius R and cross section 
7rr^, subjected to the condition r ^ R, so that the motion 
is essentially one-dimensional. The condensate is driven 
by a time-periodic potential, reading 



VExtix, t) = [Vi cos(x) -t- V2 cos(2.T + 4>)] cos [uj{t 



(1) 

where Vi and V2 denote the potential depths of the two 
lattice harmonics, </) their relative phase, and lo the angu- 
lar frequency of the oscillations of the perturbation. For 
^Li and a trap of radius R = 2.5 /im our used unit of 
time is mR^/h = 6.93 x 10~^ s. This potential can be 
readily reproduced in experiments [l8| . The space sym- 
metry (x — !■ —x-'r'-f, with 7 constant [l3|) is broken when 
(/) 7^ 0, TT. The driving is symmetric under t —t -\- 2to 
inversion and in particular, for to = 0, under t —t time 
symmetry. 

At zero temperature and as long as the number iVg 
of condensate particles is much larger than the number 
5N of non-condensate ones, the evolution of a BEG is 
well described by the Gross-Pitaevskii (GP) equation [l[. 
Thus, for our system the evolution of the condensate 
wave-function il}{x,t), normalized to 1, is described by 
(taking h = m = 1) 



dip{x, t) 



V-Eyit{x,t) 7/1(2;, t), 

(2) 

where x G [0, 27r) is the azimuthal angle of the torus, g — 
SNaR/r'^ the scaled strength of the nonlinear interaction 
(we consider the repulsive case, i.e., g > 0) and a is the 
s-wave scattering length for elastic atom-atom collisions. 

As initial condition we use a wave-function which is 
symmetric both in space, x, and in momentum, p, i.e. it 



2 



is of the form -0(2;, 0) — ane'"" cos(na;), where n G 
N and an,Oin S For simplicity, and without loss of 
generality we choose the initial condition to be: 



1 e'" 

.- cosix) 



(3) 



The asymptotic time-averaged current is shown as a 
function of the initial phase loIq of the driving field in 
Fig. [2] We can deduce that there is a strong asymmetry 
around ut^ = tt, in particular for the coupling strength 
g = 0.075. As a result, the directed current survives even 
after averaging over to (cf. the inset of Fig. (2]). 



where a is the relative phase between the homogeneous 
part and the last term which corresponds to a cosi- 
nusoidal modulation in the initial density of particles. 
While in the non-interacting (g = 0) case the asymptotic 
current is zero for any initial condition that is symmetric 
in X and p , we show that this is not the case if ^ 0. 

As depicted with Fig. [T] (top panel), our model 
may exhibit a non-zero asymptotic current, (p)asym = 
linit^oo 7 /g when ,g 7^ and a 7^ with 

= -i J^" il;*{x,t')-^ip{x,t')dx. The asymptotic 
time-averaged momentum (p)asym ff' Fig. [1] (bottom 
panel), clearly shows that (i) there is a finite threshold 
value g* such that (p)jjsyni 7^ when g > g* and (ii) there 
exists an optimal value ^opt that maximizes the current. 
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FIG. 1: (color online) Upper panel: Momentum expectation 
(p) vs. time t for a — (dark black curve) and a — tt/2 (light 
red curve). Used parameter values are: g = 0.2, Vi = V2 = 2, 
(f) = it/2, uj = 10, to = 0. Lower panel: Asymptotic time- 
average momentum {p)^^y^ versus the interaction strength g 
for a — n/2. Data are obtained from the GP equation ([2]) 
(black squares) or from the TMM-Ansatz ([5]) (red circles). 
In the first case g* ~ 0.065 and gopt ~ 0.15, in the latter 
g* ^ 0.070 and popt ~ 0.15. 
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FIG. 2: (Color online) Asymptotic time-averaged momentum 
(p)asym versus uitt) for g — 0.05 (filled black squares), g — 
0.075 (empty red circles), g — 0.1 (filled blue triangles), and 
g = 0.2 (pink asterisks). Inset: asymptotic current averaged 
over to, (p)asyi-n, as a function of the interaction strength g. 
Parameter values: Vi = 1^2 = 2, = tt/2, uj = 10, a = tt/2. 



The asymptotic time-averaged momentum is depicted 
as a function of the relative phase a in Fig. [3] (top panel). 
We note that for certain intervals of the values of a it is 
possible to obtain a non-zero current, and these regions 
widen as g increases. There is a symmetry for a lir—a 
and hence the average over all a yields zero current. 

This a —t (27r — a) symmetry can be reasoned as fol- 
lows: With io = 0, let t T = — t and perform complex 
conjugation in the GP equation ([2]), to obtain 



d'il)(x, t) 
dr 



2dx^ 



'5|-0(x,r)p -I- VEKt{x,T) 



■4j{x,t) 



(4) 

where iI){x,t) = ip*{x, —t). Therefore, if 'ip{x,t) is a 
solution of the GP equation, so is ?/'(x,t). It follows 
that (p(t)) of the wave function t) is such that 
{pi''')) — '~{p{'t))y where {p{t)) is the momentum expec- 
tation of i(j{x, t). For the initial condition ([3]), ipix, 0) — > 
ip*{x, 0) = ip{x, 0) when a (27r — a). We thus conclude 
that (p)asym(a) = - (_p) asym (27r - a) . Clearly, if -0(0) € M 



a 



0,7r in Eq. jS])] then (p). 



= 0. While in the non- 
interacting case the current is asymptotically vanishing 
for any initial condition, complex or real [l^ , in the ,9 7^ 
case a directed current can arise when a ^ 0, tt. 

In Fig. [3] (bottom panel) we show that a directed cur- 
rent can be generated only when the potential is spatially 
asymmetric, that is, 7^ 0, tt. Note that non-zero current 
is obtained and remains stable in two 0-regions around 
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FIG. 3: (Color online) Top panel: Asymptotic time-averaged 
momentum {p)asym versus a for g = 0.05 (full black squares), 
g = 0.075 (empty red circles), g = 0.15 (full blue triangles) 
and g — 0.2 (pink crosses). The other parameters values are 
Vi — V2 = 2, (f) — 7r/2, u = 10 and to = 0. Bottom panel: 
Asymptotic time-averaged momentum {p)asym versus (p for 
Vi = V2 = 0.2 (full black squares), Vi = V2 = 0.5 (empty red 
circles). Vi — V2 ~ 1 (full blue triangles) and Vi = V2 ~ 2 
(pink crosses). Other parameters values are g = 0.2, a = 7r/2, 
Lj — 10 and to = 0. 
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FIG. 4: Population of the three levels, \Af (black), |P|^(light 
gray), and \M\^ (gray) which carry momentum 0,-1-1, and -1, 
respectively, for a = (top) and n/2 (bottom), a.t g = 0.2, 
Vi = V2 = 2,(l}^ -k/2, uj = 10, to = 0. 



the maximum symmetry violation values 4> = 7r/2, 37r/2. 

From our numerical studies we notice that the only 
states that are significantly excited during the dynamical 
evolution of the system assume the momenta 0, -1-1 or 
— 1. We hence attempt to reproduce our present results, 
at least on a qualitative level, with the use of a mean-field 
three- mode model (TMM). We use the Ansatz 



1 



[A{t) + P{t)e' ^ + M{t)e- 



(5) 



where A{t), P{t) and M{t) are time-dependent, complex 
coeflacient such that |Ap + |Pp + = 1. Using ^ we 
find the effective mean-field Hamiltonian of the TMM: 



1 



(6) 



together with the equations of motion 



iA 



|p ((|^|2 -I- 2|P|2 -f 2\M\^)A + 2A*PM) 



+ f{P + M)cos{Lot), 
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iM 
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^ ((|P|2 -f- 2\A\^ + 2|A/|2)P + AHr 



-fi^ + cos(wt), 



g 

2tt 



{{\M\ 



MA- 



-2\A\^ + 2\P\^)M + A^P*) 
iP) cos{tot). 



(7) 

This TMM qualitatively reproduces the behavior of 
the asymptotic time-avcragcd current, see Fig. [T] (bottom 
panel). In Fig. 0] we depict the evolution of the popula- 
tions and |Mp of the three modes. For a = 
(a — 7r/2) the time-averaged current is zero (non-zero). 

The obtained results may appear surprising if one ob- 
serves that the full second-quantized quantum many- 
body operator is linear, obeying time-reversal symmetry, 
and symmetry arguments necessarily imply a vanishing 
asymptotic current [l^ . Actually there occurs no con- 
tradiction because the two limits t 00 and N 00 
do not commute: a directed current is obtained only in 
the mean-field limit in which we let first N 00 and 
then t —^ 00. The issue for a physical system in which 
both the number of particles and the time of the exper- 
imental run are finite constitutes the vindication of the 
regime of validity of the used nonlinear GP equation. For 
this purpose, we next study the second quantized many- 
body problem within the three-mode approximation, be- 
ing also accessible to numerical investigation. 

The second-quantized quantum many-body Hamilto- 
nian reads: 



27r 

Hes = I dx 

10 



where, within the three-mode approximation, 



1 



(at+p^e"+mte-"). 



(8) 
(9) 



a\ pt, and rh^ denote the Boson creation operators for 
the states with momenta 0, -1-1, and —1. These op- 
erators satisfy the commutation relations [a,a^] = 1, 
= 1, [m,mt] = 1 and operators corresponding 
to different modes commute. The particle conservation 
reads (a^a -I- p^p + m^rh) = N , where N is the total 
number of particles. As initial condition, we consider 



4 



atoms in the state 



V2 



+ 



N 



|0). 



(10) 



We studied numerically the evolution of the TV-body 
system governed by the second-quantized three-mode 
Hamiltonian ([8|). The normalized population inbalance 
between levels with momenta -1-1 and —1 is depicted in 
Fig. [51 The TMM holds up to increasingly longer times 
t* as N increases. We estimate t*{N) as the time instant 
at which the A'^-body population inbalance deviates by 
more than 1% from the mean-field treatment value: Our 
numerical data are well fitted by a logarithmic depen- 
dence, t* oc In N, cf. inset of Fig. [S] 

To check the stability of the condensate, we have diago- 
nalized the single-particle reduced density matrix, whose 
matrix elements read (x'^y), with x, y being a,p or rh. Ac- 
cording to [31, we have a BEG if this matrix possesses 
one eigenvalue of order N. Defining t as the time for 
which the condensed state is highly populated, we have 
found numerically that i oc Af^ "*"^. This means that for 
large and times t obeying t > t > t* the mean-field ap- 
proach no longer provides an accurate description of the 
many-body system despite the presence of a condensate. 
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FIG. 5: (Color online) Population inbalance divided by A*' 
between levels with p — +1 and p — —1 vs. time t. From 
bottom to top: N = 10, 80, and mean-field TMM, for g = 0.2, 
Vi = V2 = 2, 4> = n/2, tj = 10, to = 0, a = n/2. Inset: t* 
versus N (squares) and numerical fit t* = A + BlnN, with 
A « 73 and B « 54. 

In conclusion, we have shown that in BEG the interac- 
tion plays a prominent role in inducing long lasting cur- 
rents, when otherwise impossible. We note that BECs 
in toroidal traps have recently been realized in differ- 
ent experiments 0, [2l|, [23, [3. The initial condition 
(III) can be implemented by properly exciting the ground- 
state of a condensate in a toroidal trap, and the phase 
a can be tuned by letting the excited condensate evolve 
over a given time span. The interaction strength can be 



changed by tuning an external magnetic field [24|. Ex- 
trapolating the dependence given in Fig. [5] of t* on the 
number A^ of condensate particles, we obtain t* « 800 (in 
an in situ experiment this corresponds a time scale 0.1 
s) for N 5 X 10^ particles [l^]. Changing the number 
of condensed particles [2^ opens the possibility to ex- 
plore in situ the range of validity of the GP description 
via the measurement of our predicted, intriguing directed 
transport. 
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